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Abstract. Very recently a new inversion method has been developped to analyze 
the intergalactic medium seen in absorption in quasar spectra (the so-called Lyman- 
q forest). This method is applied to recover the temperature of the gas and the 
^ ■ underlying density field. Using constraints from the Lyman-/3 forest, it is possible 

CN ! to recover this field up to over-densities 8 ~ 10. By inverting the H I and C IV 

absorptions together it has been shown that the C iv/H I ratio varies through the 
profile of strong lines, beeing larger in the wings. The method can be applied to 
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f*^ ■ reconstruct the 3D density field from multiple lines of sight and is shown to give 

good results up to mean separations of 3 arcmin. Results from a survey of QSO 
pairs performed with HST/ST1S and VLT/UVES-FORS are summarized. 
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1 Introduction 
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-j-j ■ The numerous absorption lines seen in the spectra of distant quasars (the so- 

^ \ called Lyman-a forest) reveal the intergalactic medium (IGM) up to redshifts 

larger than 5. It is believed that the space distribution of the gas traces the 
potential wells of the dark matter. Indeed, recent numerical A-body simula- 
■ tions have been successful at reproducing the observed characteristics of the 

H ; Lyman-a forest (e.g. || f| ) . The IGM is therefore seen as a smooth all- 

pervading medium which can be used to study the spatial distribution of the 
mass on scales larger than the Jeans' length. This idea is reinforced by ob- 
servations of multiple lines of sight. It is observed that the Lyman-a forest is 
fairly homogeneous on scale smaller than 100 kpc (e.g. jn^) and highly cor- 
related on scale up to one megaparsec (e.g. The number of suitable 
multiple lines of sight is small however and the sample need to be significantly 
enlarged before any firm conclusion can be drawn (see Section 3). 
Very recently, new methods have been implemented to recover the real den- 
sity space distribution of the IGM by inversion of the Lyman-a forest |J [|J . 
It is possible to derive the physical state of the gas (temperature, density, 
metallicity) by inversion of high S/N ratio and high resolution spectra taken 
with UVES (see Section 2). The method has also been tested to recover the 
3D topology of the large scale structures when applied to the inversion of a 
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Figure 1: Part of the high resolution and high S/N spectrum of the QSO HE 
1122-1628. The normalised flux is fitted (upper panel, the solid line is the data 
and the dashed line the fit model) and the corresponding density field is re- 
covered (lower panel) . The broadening of the recovered density peaks depends 
on the temperature of the IGM. High temperature yields narrow profiles (solid 
line). Lower temperature yields smoother profiles (dashed line). 



network of adjacent lines of sight (see Section 3). 



2 Inversion of the Lyman-a forest 

We have developped a new inversion method to recover the one dimensional 
density field from the fit of normalised QSO spectraQ. An example of this 
reconstruction is given in Fig. [I] which shows part of the spectrum of HE 1122- 
1628 (upper panel) and the reconstructed density (lower panel). The fit is 
overplotted in the upper panel as a dashed line. 

To illustrate the dependence of the results on the gas temperature, the density 
field is reconstructed for different assumed constant temperatures. High tem- 
peratures yield narrow reconstructed density profiles (solid line in the lower 
panel of Fig. Q); lower temperatures yield smoother profiles (dashed line). 
Therefore, although some degeneracy exists (see M), one can constrain the 
temperature at mean density of the IGM and its cosmological evolution. 
Using the high resolution, high S/N spectrum of QSO HE1122-1628 obtained 
during Science Verification of UVES we derive a low temperature of 10 4 K at 
z ~ 2.2. A more detained description of this procedure and other applications 
of the method may be found in [|. In particular, using additional constraints 



Figure 2: Dark matter distribution in the CDM simulation, at z=2. Darker 
regions are denser (logarithmic scale) . A grid of lines of sight through the box is 
defined along the direction shown by the arrows. The inversion method is then 
applied to the synthetic spectra and results are compared to the simulation 
(see Fig. 3). 

from the Lyman-/3 forest, it is possible to recover the density field up to over- 
densities 6 ~ 10. In addition by inverting the H I and C iv absorptions together 
it has been shown that the C iv/H I ratio varies through the profile of strong 
lines, being larger in the wings ||. 

3 Large scale structures 

We have applied the method to recover the 3D spatial distribution of the IGM 
using multiple lines of sight [||. It must be emphasized that the spectral res- 
olution along each line of sight should be of the order or slightly higher than 
the corresponding spatial resolution defined as the mean separation between 
lines of sight. This means that intermediate-resolution of the kind achieved 
by VLT/FORS is enough for the corresponding observational programme pro- 
vided the blue sensitivity is large enough... of course ! 

The method has first been tested on simulations. For this we use iV-body 

SCDM simulations performed at the Institut d'Astrophysique de Paris. 

We extract a grid of lines of sight through the simulation box (see Fig. ||) 



Slice of width 1 Mpc 



Figure 3: Left panel: Spatial distribution of the DM in a slice of 1 Mpc taken 
through the simulation (Fig. 2). Other panels : reconstructed distribution 
using a grid of lines of sight separated by, respectively, 1.2, 2 and 3 arcmin. 
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Figure 4: Transverse correlation coefficient x plotted as a function of angular 
separation (panel a). Squares are the observed numbers from FORS spectra of 
QSO pairs. The solid curve is the average value from hundred pairs taken in 
the SCDM simulation, for each separation. The pair Q 0103-294A,B exhibits 
a strong correlation. This is apparent on the spectra shown in the bottom 
panel. 



and obtain the corresponding absorption spectra. The inversion method is 
then used to recover the initial 3D density field from the ID information along 
these lines of sight. Results are shown in Fig. ^ for a slice of width 1 Mpc (for 
illustrative purpose, we use 1 Mpc ~ Ih arcmin; h = H o /(100 km/s/Mpc)) 
taken through the simulation box. The left panel shows the spatial distribution 
of the dark matter as derived from the A-body simulation. The other panels 
show the reconstructed distributions using a grid of lines of sight separated by, 
respectively, 1.2, 2 and 3 arcmin. As expected, structures are better recovered 
for smaller separations. However, even with 3 arcmin separation, it is still 
possible to recover the topology of the large scale structures. 
To reach this angular separation, a density of QSOs larger than 100 per square 
degree is needed. This implies to observe QSOs of magnitude of the order of 
21.5 or fainter (see (l3)) and to chose the best field (to be found) where the 
QSO density is highest. Large surveys of QSOs are underway and such a field 
will probably be available very soon. 

In the mean time, we have started a study of the correlation between pairs of 
QSOs. From the transverse correlation it is possible to derive 3D information 
on the density field ||(7]]. Here we test the transverse correlation with a 
convenient single number, 



where J-q and T\ are the two normalised spectra of a pair of quasars separated 
by a transverse distance r±. The value of x gives a quantitative estimate of 
the correlation, ranging from if the spectra are uncorrelated to 1 for identical 
spectra. We derive x versus r± in the simulation (solid curve in Fig. g). This 
curve is model dependent mainly trough the relation between angular and 
comoving distances, but also through the size of the structures, etc... This 
method yields therefore an additional cosmological test. About twenty pairs 
are needed for each separation bin to discriminate between different cosmologic 
models @. 

We measure x on our recent FORS observations of 10 QSO pairs at z ~ 2.2 
(filled squares on Fig. 4). The observational points are consistent with the 
SCDM model. However, an interesting pair shows a much stronger correlation. 
This strong correlation is apparent from the spectra (bottom panel in Fig. 4) 
and has no explanation yet. It is also important to realise that the correlation 
is still significant at 4 arcmin. 

At smaller redshift, z ~ 1, new HST data (l5[ [l| show that a correlation signal 
is present up to 3 arcmin (1.5ft.J Mpc physical). However the coefficient x 
is smaller than at high redshift flagging expected cosmological evolution. The 
later is probably mostly due to evolution in the strength of the absorption 
lines. Comparison are difficult however because of small statistics but also 
and mostly because of comparatively smaller sensitivity of HST observations. 
Future observations with new instruments will overcome these difficulties. 
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Figure 5: Number of matches between z ~ 1 absorptions seen along the two 
lines of sight of quasar pairs with separations between 2 and 3 arcmin (~ 
1 to 1.5 h^Q Mpc physical) against the velocity difference between the two 
absorptions. The dashed line corresponds to the expected numbers for ran- 
domly distributed lines. It is apparent that there is an excess of matches for 
AV < 350 km s" 1 . 
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